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An overview of the theory of the upper critical field in dirty two-gap superconductors, with a 
particular emphasis on MgB2 is given. We focus here on the maximum H C 2 which may be achieved 
by increasing intraband scattering, and on the limitations imposed by weak interband scattering and 
paramagnetic effects. In particular, we discuss recent experiments which have recently demonstrated 
ten-fold increase of H C 2 in dirty carbon-doped films as compared to single crystals, so that the H C 2(0) 
parallel to the ab planes may approach the BCS paramagnetic limit, H P [T] = 1.84T C [K] ~ 60 — 70T. 
New effects produced by weak interband scattering in the two-gap Ginzburg-Landau equations and 
H C 2(T) in ultrathin MgB2 films are addressed. 

PACS numbers: PACS numbers: 74.20. De, 74. 20. Hi, 74.60.-w 



INTRODUCTION 



It is now well established that superconductivity in 
MgB 2 with the unexpectedly high critical temperature 
T c w 40K p]], is due to strong electron-phonon inter- 
action with in-plane boron vibration modes. Extensive 
ab-initio calculations 0, 0, 0] , along with many experi- 
mental evidences from STM, point contact, and Raman 
spectroscopy, heat capacity, magnetization and rf mea- 
surements 0, [1] unambiguously indicate that MgB 2 ex- 
hibits two-gap s-wave superconductivity 0, @|. MgB 2 
has two distinct superconducting gaps: the main gap 
A CT (0) w 7.2mV, which resides on the 2D cylindrical parts 
of the Fermi surface formed by in-plane a antibonding 
p xy orbitals of B, and the smaller gap A 7r (0) « 2.3mV 
on the 3D tubular part of the Fermi surface formed by 
out-of-plane n bonding and antibonding p z orbitals of B. 

The discovery of MgB2 has renewed interest in new 
effects of two-gap superconductivity, motivating differ- 
ent groups to take closer looks at other known materials, 
such as YNi 2 B 2 C and LuNi 2 B 2 C borocarbides @| Nb^Sn 
0, or NbSe 2 [IH, heavy-fermion and organic [13j 
superconductors, for which evidences of the two gap be- 
havior have been reported. However, several features of 
MgB 2 set it apart from other two-gap superconductors. 
Not only does MgB 2 have the highest T c among all non- 
cuprate superconductors, it also has two coexisting order 
parameters ^ a = A CT exp(i#i) and = A Tr exp(i0 2 ), 
which are weakly coupled. The latter is due to the fact 
that the a and it bands are formed by two orthogonal sets 
of in-plane and out-of-plane atomic orbitals of boron, so 
all overlap integrals, which determine matrix elements of 
interband coupling and interband impurity scattering are 
strongly reduced [14|. This feature can result in new ef- 
fects, which are very important both for the physics and 
applications of MgB 2 . Indeed, two weakly coupled gaps 
result in intrinsic Josephson effect, which can manifest 
itself in low-en erg y interband Josephson plasmons (the 
Legget mode) [15J with frequencies smaller than A % /h. 
Moreover, strong static electric fields and currents can 



decouple the bands due to formation of interband tex- 
tures of 2-7T planar phase slips in the phase difference 
6{x) = 9i — 02 flil 17 1 well below the global depair- 



ing current. In turn, the weakness of interband impurity 
scattering makes it possible to radically increase the up- 
per critical field H c2 by selective alloying of Mg and B 
sites with nonmagnetic impurities. 

Despite the comparatively high T c , the upper critical 
field of MgB 2 single crystals is rather low and anisotropic 
with H^(0) ~ 3 - 5T and H l c l 2 (0) ~ 15 - 19T of @, i], 
where the indices _L and 1 1 correspond to the magnetic 
field H perpendicular and parallel to the ab plane, re- 
spectively. Since these H c2 values are significantly lower 
than H c2 (0) ~ 30T for Nb 3 Sn 0, EH, there had been 
initial scepticism about using MgB 2 as a high-field su- 
perconductor, until several groups undertook the well- 
established procedure of H c2 enhancement by alloying 
MgB 2 with nonmagnetic impurities. The results of high- 
field measurements on dirty MgB 2 films and bulk sam- 
ples has shown up to ten-fold increase of H^> as com- 
pared to single crystals 0, ED, 0, 0, 0, 0,111, 0, IS 
0, 0, HH , particularly in carbon-doped thin films 28| 
made by hybrid physico-chemical vapor deposition [32j. 
This unexpectedly strong enhancement of H c2 (T) results 
from its anomalous upward curvature, rather different 
from that of H C 2(T) for one-gap dirty superconductors 



33j, 13J, m, \m- As shown in Fi g- 1= H c2 of MgB 2 C- 



doped films has already surpassed H c2 of NbsSn, which 
could make cheap and ductile MgB 2 an attractive mate- 
rial for high field applications 37 1 . 

This radical enhancement of H c2 shown in Fig. 1 is in- 
deed assisted by the features of two-gap superconductiv- 
ity in MgB 2 . Fig. 2 give s another example of H C 2(T) for 
a fiber-textured film [251 ] . which exhibits an upward cur- 
vature of H c2 (T) for H\\c, This behavior of H c2 (T) and 
the anomalous temperature-dependent anisotropy ratio 
r(T) = H l J 2 {T)/H^(T) are different from that of the 
one-gap theory in which the H c2 (T) has a downward 
curvature, while the slope H' c2 = dH c2 /dT at T c is pro- 
portional to the normal state residual resistivity p n , and 
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FIG. 1: H c2 {T) for carbon-doped MgB 2 films 28] in compar- 
ison with NbTi and Nb3Sn. The red and blue lines show fits 
from Eq. (QD with g = 0.045 




FIG. 2: H c2 (T) of a fiber-textured MgB 2 film g|| both par- 
allel (triangles) and parallel (squares) to the ab planes. The 
solid lines show calculations from Eq. (|f 9p with g — 0.065, 



D,r « Di ab) for H\\c and = 0.19{D { f> D, 



(c) j-,(ab),l/2 



for H _L c. 



ff c2 (0) = 0.69T c H' c2 [II, Bill, El . However, the behav- 
ior of H c2 {T) in MgB 2 can be explained by the two-gap 
theory in the dirty limit based on either Usadel equations 
H, 33] or Eliashberg equations [t| H(| . 

The behavior of H C 2(T) can be qualitatively under- 
stood using a simple bilayer model shown in Fig. 3, 
which captures the physics of two-gap superconductiv- 
ity in MgB 2 , and suggests ways by which H c2 can be 
further increased. Indeed, MgB 2 can be mapped onto a 
bilayer in which two thin films corresponding to a and 
7r bands are separated by a Josephson contact, which 
models the interband coupling. The global H c2 (T) of the 
such weakly-coupled bilayer is mostly determined by the 
film with the highest H c2 , even if Tc^ and are very 
different. For example, if the n film is much dirtier than 
the cr film then dominates at higher T, but at lower 
temperatures the 7r film takes over, resulting in the up- 




FIG. 3: The mechanism of the upward curvature of H C 2(T) 
illustrated by the bilayer toy model shown in the inset. The 
dashed curves show H C 2(T) calculated for a and n films in 
the one-gap dirty limit with the BCS coupling constants A CT = 
0.81, = 0.285, and D n = 0.LD CT . The solid curve shows 
H C 2(T) calculated from Eq. (|26[1 of the two-gap dirty limit 
theory for the BCS matrix constants from Ref. 



ward curvature of H c2 (T) . If the a film is dirtier, the 7r 
film only results in a slight shift of the H c2 curve and a 
reduction of the slope H' c2 near T c . 

The bilayer model also clarifies the anomalous angular 
dependence of H c2 (a,T) for H inclined by the angle a 
with respect to the c-axis (parallel to the film normal in 



Fig. 3) |41| . In this case both (a, T) and Hg ' (a, T) 
depend on a according to the temperature-independent 
one-gap scaling H c2 {a) = H c2 {0, T)j\J cos 2 a + esin 2 a 
(42I |43||. but with very different effective mass ratios e = 
Tn a b/m c for each film. Because the a band is much more 
anisotropic than the 7r band, e CT <C 1, and ~ 1 [44l . 
EHj . the one-gap angular scaling for the global H c2 {a,T) 
breaks down. For example, in the case shown in Fig. 3, 
H C 2(T) is anisotropic at higher T, but at lower T, the 
nearly isotropic 7r band reduces the overall anisotropy of 
H c2 , so the ratio T(T) = H l J 2 (T) / H^(T) decreases as 
T decreases. This is characteristic of many dirty MgB 2 
films like the one shown in Fig. 2, for which the 7r band 
is typically much dirtier than the a band. By contrast, 
in clean MgB 2 single crystals T(T) increases from ~ 2 — 3 



rW 



6 at t< t c M,\M,\M,\M,M,M,^2 



near T r to ~ 5 - 

5^| . This behavior was explained by two-gap 



5J, |5J, |55|, 

effects in the clean limit 

Fig. 3 suggests that H c2 (T) of MgB 2 can be signifi- 
cantly increased at low T by making the tt band much 
dirtier than the main a band. This could be done by dis- 
ordering the Mg sublattice, thus disrupting the p z boron 
out-of-plane orbitals, which form the ir band. Achiev- 
ing high H c2 requires that both a and 7r bands are in 
the dirty limit. Yet, making the 7r band much dirtier 
than the a band provides a "free boost" in H c2 without 
too much penalty in T c suppression due to pairbreak- 
ing interband scattering or band depletion due to doping 
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|59i . l60(. In fact, the interband scattering is weak for 
the same reason that ^ a and ^ n are weakly coupled, 
which may enable alloying MgB 2 with more impurities 
to achieve higher H c2 . Systematic incorporation of im- 
purities in MgB2 has not been yet achieved because the 
complex substitutional chemistry of MgB 2 is still poorly 
understood 61, 62, 6^ 64 1. Several groups have reported 
a significant increase in H C 2 by irradiation with protons 
[65j |. neutrons 6(| 6?J or heavy ions .66$, but so far the 
carbon impurities have been the most effective to provide 
the huge if c2 enhancement shown in Figs. 1 and 2. The 
effect of carbon on different superconducting properties 
can be rather complex 6^, 7^, 71 1 and still far from being 
fully understood. Yet given the indisputable benefits of 
carbon alloying, one can pose the basic question: how far 
can H C 2 be further increased? 

The bilayer model suggests that H c2 increases if intra- 
band scattering is enhanced. However, because intraband 
impurity scattering causes an admixture of pairbreaking 
interband scattering, the first question is to what ex- 
tent weak interband scattering in MgB2 can limit -ff c2 . 
Another important question is how far is the observed 
H C 2 from the paramagnetic limit H p . In the BCS theory 
Hp is defined by the condition: ^bH p cs = A/V2, or 
H^ CS [T] = 1.861; [if] [7^, where [i B is the Bohr mag- 
neton. For T c — 35-fC, this yields H p = 65T, not that 
far from the zero-field h]} 2 (0) in Figs. 1 and 2. However, 
the BCS model underestimates H p , which is significantly 
enhanced by strong electron-phonon coupling [73j: 



(1) 



Hp ~ (1 + \ ep )H p c , 

where X ep is the electron-phonon constant. Taking X ep ~ 
1 for the a band [2j, |3j, we obtain H p ~ 130T, so there 
still a large room for increasing H c2 by optimizing the 
intra and interband impurity scattering. For instance, 
increasing H' c2 to a rather common for many high field 
superconductors value of 2T/ K (much lower than H' c2 ~ 
5 - 14T/K for PbMo 6 S 8 [z|) could drive H c2 of MgB 2 
with T c ~ 35A" above 70T. In the following we give a brief 
overview of recent results in the theory of dirty two-gap 
superconductors focusing on new effects brought by weak 
interband scattering and paramagnetic effects. The main 
conclusion is that, although interband scattering in MgB 2 
is indeed weak, it cannot be neglected in calculations 
of -ff c2 (T). We will also address the crossover from the 
orbitally-limited to the paramagnetically limited iJ c2 in 
a two-gap superconductor. 



Although the a band is anisotropic, MgB 2 is not a layered 
material 75, 76| , so the continuum BCS theory is applica- 
ble because the c-axis coherence length £ c is much longer 
than the spacing between the boron planes ~ 3. 5 A. In- 
deed, even for H^(0) = 40T and H l c [(0) = 60T in Fig. 1, 
the anisotropic Ginzburg-Landau (GL) theory [42| gives 
£ c = (foH^ir) 1 / 2 jH\ k. 191. Strong coupling in 
MgB 2 should be described by the Eliashberg equations 
[4CJ , but we consider here manifestations of intra and in- 
terband scattering and paramagnetic effects in _ff c2 using 
the more transparent two-gap Usadel equations [38| 



ufi — — [3in Q n /i - /iV a V/33i] 

= ^i9i + 712(51/2 ~ 32/1) 



w/ 2 



D, 



a/3 



[92IWW2-/2V 



*232 +72l(92/l -9l/2), 



(2) 



(3) 



Here the Usadel Green's functions f m (r,u>) and g m (r,u>) 
in the m-th band depend on r and the Matsubara fre- 
quency ui — 7rT(2n +1), are the intraband diffusiv- 
ities due to nonmagnetic impurity scattering, 2 r y mm i are 
the interband scattering rates, II = V + 2iriA/(f>o, A is 
the vector potential, and 4>q is the flux quantum. Eqs. 
<j2j) and Q are supplemented by the equations for the 
order parameters ^ m = A m exp(i<p m ), 



* m = 2ttT X mm >f m < (r, uj), 



(4) 



o;>0 m 



normalization condition |/ rn | 2 
current density 



= 1, and the super- 



J a = -2neTImJ2J2 NmD "nf^pfn 



(5) 



Here N m is the partial electron density of states for both 
spins in the m-th band, and a and (3 label Cartesian 
indices. Eqs. (fj| contains the matrix of the BCS cou- 
pling constants A mm/ = A^, - fi mm / , where A^, are 
electron-phonon constants, and (J, mm ' is the Coulomb 
pseudopotential. The diagonal terms An and A 22 quan- 
tify intraband pairing, and Ai 2 and A 2 i describe inter- 
band coupling. Hereafter, the following ab initio values 
\ aa w 0.81, A™ ~ 0.285, A CT7r w 0.119, and A™ w 0.09 
[77) 1 are used. There are also the symmetry relations: 



THO-GAP SUPERCONDUCTORS IN THE 
DIRTY LIMIT 

We regard MgB 2 as a dirty anisotropic superconductor 
with two sheets 1 and 2 of the Fermi surface on which 
the superconducting gaps take the values Ai and A 2 , re- 
spectively (indices 1 and 2 correspond to a and 7r bands) . 



JV1A12 = A 2 A 2 i, A l7l2 = A 2721 (6) 

where N v k. l.ZN a for MgB 2 . Solutions of Eqs. ©-© 
minimize the following free energy J Fd 3 r [17J: 

F = \ E N ^AKL> +Fi+F 2 + F l (7) 

mm' 
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Here F\ and Fi are intraband contributions, 

F m = 2irT N m [(uj(l - g m ) - 

UJ>0 

Re{f* m A m ) + D£f [II a / m Ii;/* + V a g m W p g m ]/A 
and Fj is due to interband scattering [78j |: 

F - 2^qT - 5l52 - Re(/ 1 */ 2 )], 



(8) 



(9) 



w>0 



where 2q = iVi7i2 + -^2721- The Usadel equations result 
from 8F/8}*, = 0, dF/d^* m = 0, and J = -cSF/6A. 
Taking f m = sina m and g m = cosa m , we obtain 

wsinai + 712 sin(ai — 02) = Ai cosai, (10) 
wsina2 + 721 sm(a 2 — a i) = A2 cosa^. (11) 

These coupled equations along with Eq. (U|) define the 
two-gap uniform states for J = 0. 



CRITICAL TEMPERATURE 

Eqs. © and give the well-known results for T c 
in two-gap superconductors @, S, ZH 8(j . For negligible 
interband scattering, substitution of fi — Ai/u and = 
A 2 /lu into Eq. (4} yields: 



T c0 = 1.14frwDexp[-(A+ - A )/2w], 



(12) 



where A± = An ± A 2 2, w = AnA 22 — Ai 2 A 2 i, and A = 
(A 2 . + 4Ai2A2i) 1 ^ 2 . The interband coupling increases T c q 
as compared to noninteracting bands (A12 = A21 = 0), 
while intraband impurity scattering does not affect T c o, 
in accordance with the Anderson theorem. Solving the 
linearized Eqs. ([2]) and ((3]) with ^ mm > ^ 0, gives T c with 
the account of pairbreaking interband scattering: 



nf) = J Xo 



win t c ) hit c 



p + w In t c 

2p = A + [7-A_ - 2A21712 - 2Ai272i]/74 
17(a) = 4>(l/2 + x) -i/>(l/2), 



(13) 

(14) 
(15) 



where t c = T c /T c0 and -y± = 712 ±721, 5 = l+/2nT c0 , and 
i/j(x) is a digamma function. The dependence of T c on 
the interband scattering parameter g is shown in Fig. 4. 
As g -> 00, Eqs. (Unj) and (33) give T c -> T c0 exp(-p/io), 
and for g <C 1, we have 

7T 

T c = T c0 - -3T-[A 7+ + A-7- - 2A21712 - 2A12721] (16) 

This formula can be used to extract the interband scat- 
tering rates from the small shift of T c [8l|. However, as 
shown below, even weak interband scattering can signif- 
icantly change the behavior of H C 2(T), so it cannot be 
neglected even though g <C 1. 
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FIG. 4: Dependence of the critical temperature T c on the 
interband scattering parameter g calculated from Eq. (|13[) 
with the BCS matrix constants \ mn from Ref. XTm 



UPPER CRITICAL FIELD FOR H||c 

77 C 2 along the c-axis is the maximum eigenvalue of the 
linearized Eqs. and ([3]): 



(w ± i/j, B H)fi 



D 



Ai + (/ 2 -/i)7i2, (17) 



(w ± iii B H)f 2 - -f^h = A 2 + (h - ./ 2)721, (18) 

Here the Zeeman paramagnetic term ±/isi7, which re- 
quires summation over both spin orientations in Eq. 
(3), is included. In the gauge A y — Hx, the solu- 
tions are f m (x) = f m exp(-7ri7a; 2 /0o), and A_ m (x) = 
A m exp(— ttHx 2 / 4>q), where f m is expressed via A m from 
Eqs. (JT7J) and The solvability condition g]) of two 

linear equations for Ai and A2 gives the equation for i7 C 2 
[38| , which accounts for interband and intraband scatter- 
ing and paramagnetic effects: 



(A + XJQnt + U+) + (A - XJQnt + U-) 
+2w(lat + U + )(lnt + U-) = 0, 

where t = T/T c q, and 

Xi = [(oj- + 7_)A_ - 2A12721 - 2A 2 i7i 2 ]/fio, 
2Q± = uj + + 7+ ± f2 , 

n = [(w- + 7-) 2 + 4712721] 1/2 , 

w± = {D 1 ±D 2 )irH/4> , 



U± 



Re</>( - 



1 fi H 



2irT 



(19) 



(20) 
(21) 
(22) 
(23) 

(24) 



If interband scattering and paramagnetic effects are neg- 
ligible, Eqs. (fl"9|) - (|24"l) reduce to a simpler equation 
38l . |39j , which can be presented in the parametric form: 



lnt = -[U(h) + U{T]h)+\o/w]/2+ (25) 
[(U(h) - U( V h) - A_/ w ) 2 /4 + A 12 A 2 i/ W 2 ] 1/2 , 

77 c2 = ZfaTcth/Dt, (26) 
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where 77 = D2/D1, and the parameter h runs from to 00 
as T varies from T c to 0. For equal diffusivities, 77 = 1 , Eq. 
([2"6")> simplifies to the one-gap de-Gennes-Maki equation 
\nt + U(h) = [II SI SI]. 



Now we consider some limiting cases, which illustrate 
how H c2 depends on different parameters. Fig. 5 shows 
the evolution of H c2 as g increases for fixed D x and D 2 
and negligible paramagnetic effects, fnterband scattering 
reduces the upward curvature of H c2 (T), H c2 (0), and 
T c , while increasing the slope H' c2 at T c . Notice that 
the significant changes in the shape of H c2 {T) in Fig. 5 
occur for weak interband scattering (g <C 1), which also 
provides a finite H c2 (0) even if D 2 — > 0. For example, 
the high-field films in Fig. 1 and 2 have g ~ 0.045 and 
0.065, respectively. For g < 1, Eq. (T5]) yields the GL 
linear temperature dependence near T c : 



8 fa (T c - T) 

7T 2 (siL>l + S 2 D 2 ) 



(27) 



where T c is given by Eq. (|T6|) . si = 1 + A_/A and 
S2 = 1 — A_/Ao- Eq. I|27p is written in the linear accuracy 
in g <C 1. Higher order terms in g not only shift T c but 
also increase the slope H' c2 at T c , as evident from Fig. 5. 
For si ~ s 2 , the slope i?^ 2 is mostly determined by the 
cleanest band with the maximum diffusivity. However, 
because of weak interband coupling in MgB2, the values 
of s-y and s 2 are very different. For An = 0.81, X 22 = 
0.285, A 12 = 0.119, A 21 = 0.09 [z3, we get A_ = A n - 
A 22 = 0.525, A = (A 2 . + 4Ai 2 A 2 i) 1/2 = 0.564, thus s x = 
1 + A_/A = 1.93, s 2 = 1 - A_/A = 0.07. Thus, H' c2 
is mostly determined by D\ of the a band. Yet, if the a 
band is so dirty that Di/D 2 < s 2 /si ~ 0.04, the slope 
H' c2 is determined by the much cleaner n band. 

At low T both the Zeeman and interband scattering 
terms in Eq. (fl"9|) can be essential. Eq. (fl~9)) reduces to 
the following equation for H c2 (0): 



(A 



= w In 



,,2 Tj2 



n 2 + 



+ (A - A*) In ■ 



in - 



(28) 



where fJ-sHp = ttT c q/2j is the field of paramagnetic in- 
stability of the superconducting state, and In 7 = 0.577. 
We first consider the limit g — > 0, which defines the max- 
imum -ff C 2(0) achievable in a dirty two-gap superconduc- 
tor with no T c suppresion. In this case fi+ = nDiH / 'fa 
and f2_ = irD 2 H/fa, so for T <C T c , paramagnetic effects 
just renormalize intraband diffusivities in Eq. 



D m ^D m = ^Dl + Dl (29) 
where Dq = /igfa/n is the quantum diffusivity 

D = h/2m, (30) 
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FIG. 5: Effect of interband scattering and the diffusivity ratio 
on the evolution of H C 2(T). Upper panel shows Hc.2{T,rj) for 
the fixed g = 0.03 and different g = D 2 /D x : 0; 0.05; 0.1; 0.5 
(from top to bottom curves). Lower panel shows H C 2(T, g) for 
the fixed D2/-D1 = 0.05 and different g = 0.01; 0.05; 0.1; 0.5 
(from top to bottom curves). 



and m is the bare electron mass. Eq. (130(1 follows from 
the basic diffusion relation I 2 = Dot, and the energy un- 
certainty principle h 2 /2ml 2 — h/t for a particle confined 
in a region of length I. For g — 0, Eq. (f2"8"|) yields 



H c2 (0) = 



27 VlhDi 



exp({), 




(31) 



— . (32) 



Ao 

w 



If D < D m , Eqs. dnj-© reduce to the result of Rcf. 
(38|, and for the symmetric case, Di = D 2 , Eqs. (|3T1) - 
([52")) give the one-band result ^2(0) = faT c /2jD^^. 
However for D\ ^ D 2 , H c2 (0) can be much higher than 
H c2 (Q) — 0.69H' c2 T c . Indeed, if the effective diffusivities, 
£>i and D 2 are very different, Eqs. -([221) yield 



H c2 (0) = 
H c2 {0) 



t>oT c _ 



2jD 2 



(\-+\q)/2w 



IJ _0^ e -(Xo-\-)/2w 



2 1 D 1 



D 2 < D x e 



(33) 
(34) 



Thus, ir c 2(0) is determined by the minimum effective 
diffusivity, but unlike the limit Do — > 0, H c2 (0) remains 
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FIG. 6: Crossover from the orbitally to paramagnetically 
limited H c2 (T) calculated from Eqs. (J26J) and j37j-([38]) for 
D 2 = 0.05Di and D 1 /D = 0,1,5,20 from top to bottom 
curves, respectively. 



finite even for D\ — > or D2 — > 0. In fact, if both 
D\ -C -Do an( l -D2 ^ we return to the symmetric 
case D\ = D2, for which Eqs. (f3TT) -(f32 |) yield the result 
of the one-gap dirty limit theory [36j 



H c2 (0) ^H p = </> T c /2 7 A) = 7rr c /2 7Ms 



(35) 



For a one-band superconductor, Eq. (f35|) can also 
be written as the paramagnetic pairbreaking condition, 
fi B H p = A(0)/2, where A(0) = ttT c /7 is the zero- 
temperature gap. For two-band superconductors, the 
meaning of H p is less transparent, yet the maximum H p 
expressed via T c is given by the same Eq. (|35|) as for 
one-band superconductors. 

Finally we consider how paramagnetic effects affect the 
shape of iJ c2 (T) in the limit g — > 0. This case is described 
by Eq. ([26)1 modified as follows: 



U(h) -> Ret/' [1/2 + h(i+p)] - ip{l/2) 
U{rjh) -> ReV[l/2 + h(jrq + i)} - i>(l/2) 
H c2 = H c2 = 2<f> T c th/D , 



(36) 
(37) 
(38) 



where p = Dx/Do, and rj = D2/-D1. Fig. 6 shows how 
-Hc2 evolves from the orbitally-limited H c2 (T) with an 
upward curvature at D\ 3> D 2 to the paramagnetically- 
limited H C 2(T) of a one-gap superconductor for D\ < Do 
[72I ]. The nonmonotonic dependence of H c2 (T) in Fig. 6 
indicates the first order phase transition, similar to that 
in one-gap superconductors. 



THIN FILMS IN A PARALLEL FIELD 

H C 2 can be significantly enhanced in thin films or mul- 
tilayers, in which MgB2 layers are separated by nonsuper- 
conducting layers. It is well known that in a thin film of 
thickness d < £ in a parallel field, = 2y/SH c2 ^/d 



can be higher than the bulk H c2 = 4>o/2n£, 2 [7fl[82j]. Let 
us see how this result is generalized to two-gap super- 
conductors. For a thin film of thickness d < max^,^), 
the functions fx and / 2 are nearly constant, so integrat- 
ing Eqs. (T7|) and US} over x with d x f{±d/2) = 0, re- 
sults in two linear equations for fx and f 2 with II 2 = 
(irHd/4> ) 2 /3. Thus, we obtain the previous Eq. flTJ]) - 
([24l in which one should make the replacement 



,(/) 



{*Hd/<M 2 {D 1 ±D 2 )/& 



(39) 



We first consider the case of negligible interband scatter- 
ing and paramagnetic effects. Then Eq. (fl"9[) and (|39|) 
give the square-root temperature dependence near T c 



H. 



(/) _ A(p^iT c (T c -T) 

' TT^diSxDx + S2D2) 1 / 2 



(40) 



characteristic of thin films [82| instead of the bulk GL 
linear dependence (|27j) . From Eqs. (|3T|) and (|39|) we can 
also obtain H^P(0) for D < D TO : 



H 



(/) 

c2 



(0) = £ ( ^ 



1/2 exp(//4) 
(M2) 1 / 4 



(41) 



Next we consider the crossover to the paramagnetic limit 
in thin films at low temperatures. For neglect interband 
scattering, the expressions H%H 2 + il 2 under the loga- 
rithms in Eq. (^SJ) become fi 2 B H 2 + (ir H d / faY D 2 /36. 
Substituting here H^) ~ 4>o/£,d, we conclude that para- 
magnetic effects become essential if 



mm(D 1 ,D 2 ) < DoZ/d. 



(42) 



Thus, reducing the film thickness extends the region of 
the parameters where H C 2 is limited by the paramagnetic 
effects rather than by impurity scattering. 



ANISOTROPY OF H c i AND H c2 

For anisotropic one-gap superconductors, the angular 
dependence of the lower and the upper critical fields is 
given by [42l[43l| 



H cl (a,T) 



H cl (0,T) 
R{a) : 



H c2 (a,T) 



H c2 (0,T) 
R{a) 



(43) 



where R(a) = (cos 2 a + esin 2 a) 1 / 2 , e — m a b/m c . Here 



the anisotropy parameter T(T) = Hj^/H^ = e ' is 
independent of T for both H c \ and H c2 . By contrast, 
r 2 (T) = H l J 2 /H^ 2 for MgB 2 single crystals increases from 



c2 

- 2 - 3 at T c to 5 - 6 at T < T c , but Tx(T) = H l JjH^ 
decreases from ~ 2 — 3 to ~ 1 as T decreases [46|, |47, 

This behavior was 



48, 4£ 



51, 5/ 



54, 5E 



explained by the two-gap theory in the clean limit [57 
58, M, 0. 
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The dirty limit is more intricate in the sense that T(T) 
can either increase or decrease with T, depending on the 
diffusivity ratio D 2 /Di. However, the physics of this de- 
pendence is rather transparent and can be understood 
using the bilayer toy model as discussed in the Intro- 
duction. Indeed, for very different D\ and Z?2, both the 
angular and the temperature dependencies of H C 2(a,T) 
are controlled by cleaner band at high T and by dirtier 
band at lower T. For instance, if D2 -C D\, the high- 
T part of H C 2(a,T) is determined by the anisotropic a 
band, while the low-T part is determined by the isotropic 
7r band. In this case T(T) decreases as T decreases, as 
characteristic of dirty MgB2 films represented in Figs. 1 
and 2. If the ir band is cleaner than the a band, T(T) 
increases as T decreases, similar to single crystals. 

For the field H inclined with respect to the c-axis, the 
first Landau level eigenfunction no longer satisfies Eqs. 
([17]). fT8|) and Q. In this case f m (w,r) are to be ex- 
panded in full sets of eigenfunctions for all Landau levels, 
and H r2 becomes a root of a matrix equation M(H c2 ) = 
[HI, [33| ■ As shown in Ref. [38| , this matrix equation for 
H C 2 greatly simplifies for the moderate anisotropy char- 
acteristic of dirty MgB2 for which all formulas of the 
previous section can also be used for the inclined field as 
well by replacing D\ and D2 with the angular-dependent 
diffusivities D\{a) and D 2 {a) for both bands: 



D m (a) = [D$ 3 cos 



(44) 



In terms of the bilayer model shown in Fig. 1, Eq. (|44p 
just means that Eq. (|43p should be applied separately 
for each of the films. For g = 0, Eqs. (J57J) and (f44|) 
determine the angular dependence of H C 2(a) near T c , and 
the London penetration depth k a p is given by [38} 



A-^^^A.tanh^ 



N 2 D 



Q/3 A 2 tanh^l 



2T J 

(45) 

Eqs. (|27"|). (|4*4l . and (J60j) show that the one-gap scaling 
(|43| breaks down because the behavior of H c i(a,T) is 
mostly controlled by the cleaner band for all T, while the 
behavior of H c2 (a, T) is determined by the cleaner band 
at higher T, and by the dirtier band at lower T. Thus, 
Ti(T) and T 2 (T) for H cl and H c2 in the two-gap dirty 
limit are different. Temperature dependencies of T(T) 
were calculated in Refs. |38l.l39}. 

Eqs. (|4"4"1) and (fH)]) describe well both the temperature 
and the angular dependencies of H C 2(a, T) in dirty MgB2 
films m,El,[2i|[3i|. Eq. flS) is valid if the a band is not 
too anisotropic, and the off-diagonal elements M mn ~ 
C rn+n are negligible provided that C < 1 & Here 



c 



(ei - e 2 ) 2 sin 4 a 



{■<</ cos 2 a + e\ sin 2 a + y cos 2 a + e 2 sin 2 a] 



(46) 



ei = D[ c) /D[ ab) and e 2 = D ( 2 c) / D ( 2 ab) . For e 2 = 1, the 
parameter C( a ) < 0-45 for a rather strong anisotropy 



€\ < 0.04 and a — tt/2. For a stronger anisotropy, the 
condition £(a) <C 1 can still hold in a wide range of a, 
except a vicinity of a « i/2. In this case the calcu- 
lation of H C 2(a,T) requires a numerical solution of the 
matrix equation for H C 2 [39}. However, the Usadel the- 
ory can only be applied to dirty MgB 2 samples which, 
contrary to the assumption of Ref. [39} . usually ex- 
hibit much weaker anisotropy (r2 — 1 — 2) than single 
crystals. Perhaps, strong impurity scattering and ad- 
mixture of interband scattering reduce the anisotropy 
of D[ c) /D[ ab) ~ 0.2 - 0.3 as compared to that of the 
Fermi velocities (w 2 ) a j (f 2 b ) „ ~ 0.02 predicted by ab- 
initio calculations for single crystals 44|. The moderate 
anisotropy of D\ in dirty MgB2 makes the scaling rule 
(|4"4")) a very good approximation, as was recently con- 
firmed experimentally [30]. 



GINZBURG-LANDAU EQUATIONS 

The two-gap GL equations were obtained both for the 
dirty limit without interband scattering [38l . [85| , and for 
the clean limit [86j. Here we consider the GL dirty limit, 
focusing on new effects brought by interband scattering. 
For 7 mm ' = 0, the Usadel equations near T c yield 

f m = * m /w + D ma U 2 a */2w 2 - ¥ m \* m \ 2 /2u 3 , (47) 

where the principal axis of D a p are taken along the crys- 
talline axis. For weak interband scattering, the free en- 
ergy F = Fo + Fi contains the free energy Fo{^i, iff 2} for 
7mm' — an d the correction F^{^i, iff 2 } linear in 7 mm y. 
Here F does not have first order corrections in 7„ lm < if 
^> m satisfies the GL equations, so Fi can be calculated 
by substituting Eq. fl7|) into Eq. §§§ and expanding 

5m«l-|/m| 2 /2-|/m| 4 /8: 



F i =irqTY,[\h-f2? 



(l/il 2 



l/ 2 | 2 ) 2 /4], 



(48) 



u>0 



where q = (AT 1 7 12 + -/V2721 )/2. Combining Fi with F in 
the dirty limit for g = [381 ] . we arrive at the GL free 
energy J FdV for g <C 1: 



F = a^^ 2 + c^ll^] 2 + b 1 \y 1 \ i /2 

+ a2 |^ 2 | 2 + C 2 Q |n Q * 2 | 2 + 6 2 |^2| 4 /2 



A-I*il 2 l*2| 2 



•2&i(|*i| 2 + |* 2 | 2 )Re(*i* 2 ) 



Here the GL expansion coefficients are given by 

7T712] 



«1 



0,2 



N 1{ T 

— In ttt 

2 L T x 



2 r, -l 

— In 777- 

2 L T 2 



Clc 



f-'2o 



NiD la [ 
: N 2 D 2a [ 



7T 

16T " 
16T 



4T J 

7T721 -I 

4T ^ 

7C(3)7l2 n 
87r 2 T 2 J 

7C(3)72i - 
8tt 2 T 2 ■ 



(49) 

(50) 
(51) 
(52) 
(53) 
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h = N 1 [ 
b 2 =N 2 f 



7C(3) 3tt 7i 



16tt 2 T 2 384T 3J ' 

7C(3) 37T 721 I 



T6tt 2 T 2 384T 3J ' 
Ni 1-A12 , fl-712-1 , iVa r A 2 i , 7T721-, 
2 ^ u> 
7C(3) 



(4ttT) 



4T J 2 L w 4T J ' 
2 (Di +£> 2 )(7i2^Vi +721^2), 



6,: = 



7T 



384T 3 



(712-W1 + 721-^2), 



(54) 
(55) 
(56) 
(57) 
(58) 



where Ti = T c oexp[— (Ao — X—)/2w], and T 2 = 
T c o exp[— (Ao + A_)/2u>]. The GL equations are obtained 
by varying J FdV. I would like to point out the misprints 
with wrong signs of <2j, c\ and c 2 in Eqs. (13), (14) and 
(20) in Ref. (see also Ref. Q). 

The first two lines in Eq. (|49[) are the GL intra- 
band free energies and the term aiRe^i^) describes 
the Josephson coupling of ^1 and ^ 2 . Inter band scat- 
tering increases a\ and a 2 , and the interband coupling 
constant at- The net result is the reduction of T c deter- 
mined by the equation Aai(T c )a 2 (T c ) — a 2 , which repro- 
duces Eq. (|16p . Besides the renormalization of a m , b m 
and c m , interband scattering produces new terms, which 
describe the mixed gradient coupling and the nonlinear 
quatric interaction of ^1 and ^ 2 . Similar terms were in- 
troduced in the GL theories of heavy fermions [l2| and 
borocarbides [89j], and phenomenological models of H c2 
in MgB2 [87| • These terms result from interband scatter- 
ing, so both Ci and bi vanish in the clean limit [86| . The 
mixed gradient terms in Eq. (|49j) produce interference 
terms in the current density J = —cSF/SA.: 



r = -[(2c x A 2 + CjAiAa cos6»)Qi + 
(2c 2 Aj + c t AiA 2 cos6»)Q 2 + 
c 4 (A 2 VAi - AiVA2)sin6»]27rc/0o 



(59) 



where Q m = V9 rn + 27rA/0 o , and = 9\ — 9 2 . Here 
J is no longer the sum of independent contributions of 
two bands, because phase gradients in one band pro- 
duce currents in the other. Moreover, J acquires new 
cos 9 terms and the peculiar sin (9 interband Josephson- 
like contribution for inhomogeneous gaps. For currents 
well below the depairing limit, both bands are phase- 
locked (9 = 0), and Eq. (f5"9")) defines the London pene- 
tration depth A 2 = ccj>oQ/8ir 2 \J\: 

A = o /4^[27r( Cl A 2 + e;AiA 2 + c 2 A 2 )] 1/2 (60) 

where c\, c 2 and Cj depend on the field orientation ac- 
cording to Eq. (HH). Eq. (|47|) can be used to calculate 
H^(T) from the linearized GL equations, which give H c2 
as a solution of the quadratic equation [89( 



2i:ciH 



ai 



2irc 2 H 

4>o 



a 2 



2TTCjH 



00 



which reduces to Eq. (|27| near T c to the linear accuracy 
in 7mm'- However, GL calculations of H c2 {T) in MgB 2 
beyond the linear T c — T term [§3, [HI have a rather lim- 
ited applicability, since a\(T) and a 2 (T) change signs at 
very different temperatures T\ and T 2 . For A mn of Ref. 
[77j.7i ~ 0.9T c0 and T 2 - 0.1T c0 so higher order gradient 
terms (automatically taken into account in the Eliash- 
berg/Eilenberger/Usadel based theories) become impor- 
tant. For example, at T w T\ where a2(Ti) 3> ai(T), 
retaining the first gradient term oc c 2 requires taking into 
account a next order term ~ H 2 in the first brackets in 
Eq. (|6"Tj) . which is beyond the GL accuracy. Thus, apply- 
ing the GL theory in a wider temperature range [871 . 1881 ] 
makes it a procedure of unclear accuracy, which can re- 
sult in a spurious upward curvature in H c2 (T) not always 
present in a more consistent theory (for example, in the 
dirty limit at D\ ~ D 2 ). In addition, the anisotropy of 
Di(a) may further limit the applicability of the GL the- 
ory for H||a6, as for ci 3> c 2 higher order gradient terms 
in the ir band become important [HI . 



DISCUSSION 

The remarkable ten-fold increase of H c2 (T) in C-doped 
[H, [27|, [H, El] has brought to focus new 



MgB 2 films [25 



(61) 



and largely unexplored physics and materials science of 
two-gap superconducting alloys. Moreover, the observa- 
tions of H c2 close to the BCS paramagnetic limit poses 
the important question of how far can H c2 be further 
increased by alloying. This possibility may be naturally 
built in the band structure of MgB2 , which provides weak 
interband coupling and weak interband scattering, thus 
allowing MgB2 to be alloyed without strong suppression 
of T c . For example, for the C-doped MgB2 film shown in 
Fig. 1, p n was increased from ~ 0.4/ificm to 560/ificm, 
yet T c was only reduced down to 35K 28] . It is the weak- 
ness of interband scattering, which apparently makes it 
possible to take advantage of very dirty tt band to signif- 
icantly boost H c2 in carbon-doped films which typically 
have -Dtt ~ 0.lD a . The reasons why scattering in the 7r 
band of C-doped MgB2 films is so much stronger than 
in the a band has not been completely understood, but 
another immediate benefit for high-field magnet applica- 
tions [13] is that carbon alloying significantly reduces the 
anisotropy of H c2 down to T(T) ~ 1 — 2. 

Despite many yet unresolved issues concerning the two- 
gap superconductivity in MgB 2 alloys, H c2 of C-doped 
MgB2 has already surpassed H c2 of NbaSn (see Fig. 
1). Given the intrinsic weakness of interband scattering, 
which enables tuning MgB 2 by selective atomic substitu- 
tions on Mg and B sites, there appear to be no funda- 
mental reasons why H c2 of MgB2 alloys cannot be pushed 
further up toward the strong-coupling paramagnetic limit 
(TT|). Thus, understanding the mechanisms of intra and in- 
terband impurity scattering in carbon-doped MgB2, and 
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the competition between scattering and doping effects 
becomes an important challenge for the computational 
physics. For instance, it remains unclear why the multi- 
phased C-doped HPCVD grown films [32| exhibit higher 
H C 2 and weaker T c suppression [28j than uniform car- 
bon solid solutions 26|, 69|, 7(J. This unexpected result 
may indicate other extrinsic mechanisms of H C 2 enhance- 
ment, which are not accounted by the simple two-gap 
theory presented here. Among those may be effects of 
electron localization or strong lattice distortions in mul- 
tiphased C-doped films which can manifest themselves 
in the buckling of the Mg planes observed in the dirty 
fiber-textured MgB2 films shown in Fig. 2 [5^1. Such 
buckling may enhance scattering in the it band formed 
by out-of-plane p z boron orbitals. 

Recently significant enhancements of vortex pinn ing 
and critical current densities J c in MgB2 [90l. l9ll. 92. 93. 
[94], [95|, [9(1 has been achieved, particularly by introduc- 
ing SiC ![92j and ZrB2 [§5| nanoparticles. Given these 
promising results combined with weak current blocking 
by grain boundaries (9?]], the lack of electromagnetic 
granularity 98[ , and very slow thermally-activated flux 



creep 9{| 100l ]. it is not surprising that MgB 2 is being 



regarded as a strong contender of traditional high-field 
magnet materials like NbTi and NbaSn. Despite these 
achievements, a detailed theory of pinning in MgB2 is 
trill lacking. Such theory should take into account a 
composite structure of the vortex core, which consists of 
concentric regions of r adius £ g and where Aj(r ) and 
A x (r) are suppressed [lOl| [l02|, [Tol [Tol \}M, Em ■ For 
example, in MgB2 single crystals the larger vortex cores 
in the n band start overlapping above the "virtual up- 
per critical field" H v — o /27r£^ ~ 0.5T, causing strong 
overall suppression of A„- well below H C 2 [iMIIqI- This 
effect can reduce J c at H > H v , however both H v and 
H C 2 can be greatly increased by appropriate enhancement 
of impurity scattering in MgB2 alloys. 

Recently there has been an emerg i ng in terest in mi- 

and a possi- 



108 



crowave response of MgB2 [107 
bility of using MgB 2 in resonant cavities for particle ac- 
celerators [lid llllj . These issues require understanding 
nonlinear electrodynamic s and cur rent pairbreaking in 
two-gap superconductors 112l 1 13 1 . in particular, band 
decoupling and the formation of interband phase textures 
at strong rf currents H, 17|. 

This work was partially supported by in-house research 
program at NHMFL. NHMFL is operated under NSF 
Grant DMR-0084173 with support from state of Florida. 
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